Neutrino-nucleus interaction rates at a low-energy beta-beam facility 



"\f : 
o ■ 
o ■ 
<n : 
-«-> ■ 

o ■ 

O: 
(N : 

(N ■ 

> : 

m . 
o\ ■ 

(N ■ 

co : 
o . 
^t- ■ 
o ; 

^ : 

i 

Q . 
■ 



X 



Julien SerrearQ 

Institut fiir Theoretische Physik der Universitat Heidelberg, 
Philosophenweg 16, D-69120 Heidelberg, Germany 

Cristina VolpcQ 

Institut de Physique Nucleaire, F-91406 Orsay Cedex, France 
(Dated: February 1, 2008) 

We compute the neutrino detection rates to be expected at a low-energy beta-beam facility. We 
consider various nuclei as neutrino detectors and compare the case of a small versus large storage 
ring. 

PACS numbers: 25.30.Pt, 26.50.+X 



I. INTRODUCTION 

The pioneering experiment of R. Davis [l| has started 
the era of neutrino astronomy. Because they only have 
weak interactions with matter, neutrinos are precious 
messengers of what happens in the interior of stars, like 
our sun, or in explosive phenomena, such as Supernova 
type II explosions. Such astronomical neutrinos there- 
fore provide an important source of information for our 
understanding of the life and death of stars. Nuclei are 
commonly used as detectors in neutrino observatories as 
well as in various experiments aiming at studying intrin- 
sic neutrino properties, such as their masses and mixings. 
A precise knowledge of neutrino-nucleus cross-sections 
is needed for the interpretation of these measurements 
and/or to study the feasibility of new projects. The 
understanding of neutrino-nucleus interactions is also of 
crucial importance for various astrophysical processes. 
Timely examples include neutrino nucleosynthesis 0, 01 , 
or the nucleosynthesis of heavy elements during the so- 
called r-process 0, IE IE SB- If the latter takes place dur- 
ing the explosion of Supernovae type II, where a gigantic 
amount of energy is emitted as neutrinos of all flavors, 
final abundances depend on several nuclear properties, 
among which the interactions with neutrinos. 

According to existing simulations, the average energy 
of neutrinos emitted from core-collapse Supernovae is 
about 10 MeV for electron neutrinos and about 20 MeV 
for muon and tau neutrinos Notice however that, 
due to oscillations, electron neutrinos can become hotter 
while traversing the star 

[llHillll. The predicted spec- 
tra cover the 50 MeV region and present a tail up to about 
100 MeV 0- Reactor and solar neutrinos have typical 
energies in the 10 MeV energy range, while accelerator 
and atmospheric neutrinos cover the GeV and multi-GeV 
range. The various theoretical approaches employed to 
describe neutrino-nucleus interactions therefore involve 
nuclear as well as nucleonic degrees of freedom (for a re- 
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view, see [HQ). There are a number of open issues 
in this context. The A=2 system is the simplest case, 
for which the reaction cross sections can be estimated 
with high accuracy |15| . However, there is still an impor- 
tant quantity, namely I^i.a, related to the axial two-body 
current, which dominates the theoretical uncertainty in 
neutrino-deuteron interactions. For heavier nuclei, in the 
tens of MeV energy range, the reaction cross sections are 
dominated by collective modes, like the Gamow- Teller 
resonance or the Isobaric Analog State, which have been 
extensively studied in the past [ijj . As the neutrino im- 
pinging energy increases, transitions to states of higher 
multipolarity (such as the spin-dipole or higher forbidden 
transitions) become important [jj]]- The latter also play 
an important role in the context of core-collapse Super- 
nova physics H El El [H . Although some information 
on these states can be gathered th roug h other probes, 
such as charge- exchange reactions |l6l. muon capture 
[2fj| , or inelastic electron scattering |2l|. the experimental 
information is rather scarce. Note that the understand- 
ing of neutrino-carbon reactions with neutrinos produced 
from the decay in flight of pions is still an open issue, 
for most of the theoretical calculations over-estimate the 
experimental value 22] . So far, measurements with low- 
energy neutrinos have been performed in a few cases only, 
namely deuteron [2^ , carbon [24[ , and iron [25J . System- 
atic studies would be of great importance both for what 
concerns the interpolation from the MeV to the GeV neu- 
trino energy range and the extrapolation to neutron-rich 
nuclei, as required in the astrophysical context. 

Neutrino-nucleus interaction studies were one of the 
main physics issues of the proposed ORLAND under- 
ground neutrino facility, which was based on a conven- 
tional neutrino source (pion and muon decays) HHH. 
A smaller version of the ORLAND project is now un- 
der study 27]. At present, the MINERVA project 
includes the study of neutrino-nucleus interactions for 
neutrino energies in the GeV range. Here, we study 
the potential of a low-energy neutrino facility based on 
beta-beams, a novel method to produce neutrino beams 
|29|. This consists in boosting exotic ions which decay 
through beta-decay and produce pure, collimated and 
well-understood electron neutrino fluxes. Such a method 
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could be exploited for a future facility at CERN [29 H 30 | . 
High energy beta-beams would be fired to a gigantic 
Cherenkov detector like UNO j^H], located in an (up- 
graded) Frejus underground laboratory to study, in par- 
ticular, the possible existence of CP violation in the lep- 
tonic sector |2£j, • The discovery potential with a 
very high 7 and a longer baseline is discussed in [33l l34| . 

It has recently been proposed to use the beta-beam 
concept for the production of low-energy neutrinos [35j . 
Several laboratories will produce intense exotic beams in 
the near future and could, therefore, be possible sites for 
a low-energy beta-beam facility. These include GANIL, 
CERN, GSI, as well as the EURISOL project. Low- 
energy neutrino beams would offer an interesting oppor- 
tunity to study various neutrino properties, such as e.g. 
the neutrino magnetic moment [3q . as well as neutrino- 
nucleus interactions, of interest for nuclear physics, par- 
ticle physics and astrophysics. In the former case, one 
would exploit the ions at rest as an intense neutrino 
source, whereas, in the latter case, one would use boosted 
ions, which would be collected in a storage ring [3{|, as 
in the original high energy proposal. An important fea- 
ture of such beta-beams is that the boost factor of the 
accelerated ions can be varied, allowing one to explore 
various neutrino energy ranges. 

In this paper, we present for the first time charged- 
current neutrino-nucleus interaction rates achievable at a 
low-energy beta-beam facility. We consider two possible 
cases for the dimensions of the storage ring, for which 
we inspire ourselves of the one planned in the future GSI 
facility 1371 and the one thought in the CERN baseline 
scenario |29l l3Cj . We consider various target nuclei as 
neutrino detectors, namely deuteron, oxygen, iron and 
lead, which are commonly used in existing or planned 
experiments 1261 . Related work in the case of lead can 
be found in 1381. 



II. FORMALISM 



the emitted lepton (electron or positron), Q being the Q- 
value of the reaction, and the Fermi function F(±Z, E e ), 
which accounts for the Coulomb modification of the spec- 
trum. 

In the laboratory frame, where the boosted nucleus has 
a velocity v = (3c, the decay rate reads: 



dW 
~dT 



lab 7 47T 



(3) 



where 7 = 1/y 1 — /3 2 is the time dilation factor and 
where E v and Q = (9, ip) now denote the energy and 
solid angle of the emitted (anti-)neutrino in the labo- 
ratory (lab) frame, 9 being the angle of emission with 
respect to the beam axis. The boosted flux <^i a i,{E v ,ff) 
is given by: 



$ cm (^7[l-/3cosfl]) 
7[l-/?cos0] 



(4) 



We consider a storage ring of total length L with a 
straight sections of length D. In the stationary regime 
the mean number of ions in the storage ring is 7T5, where 
r = ^i/2/ m 2 is the lifetime of the parent nuclei and g 
is the number of injected ions per unit time. The total 
number of neutrinos emitted per unit time from a portion 
di of the decay ring is 



dN v 



dW 



dt 1T9X ^ 



X 

lab L 



(5) 



For simplicity, we consider a cylindrical detector of ra- 
dius R and depth h, aligned with one the straight sec- 
tions of the storage ring, and placed at a distance d from 
the latter. After integration over the useful decay path 
and over the volume of the detector, the total number of 
events per unit time is: 



dN m 
dt 



= grnh x 



/>oo 

/ dE v $ to t(E v )<T(E„), 
Jo 



(6) 



A. Neutrino fluxes and interaction rates 

The decay rate of a nucleus in the rest (cm) frame can 
be written as: 



dW 
~dT 



d 2 tt 

<$> cm (E v )dE v —— , 

4-7T 



(1) 



where E v and £1 denote respectively the energy and the 
solid angle of the emitted (anti-)neutrino, and where the 
neutrino flux <& cm (E v ) is given by the well-known formula 
H3: 

®cm{E v ) =bElE e y/E* - ml F(±Z, E e ) ®{E e - m e ) . 

(2) 

where the constant b = In 2/m^,ft 1 /2, with m e the elec- 
tron mass and ft 1 / 2 the ft- value. The quantities appear- 
ing in the above expression are the energy E e = Q — E^oi 



where n is the number of target nuclei per unit volume, 
a(Ev) is the relevant neutrino- nucleus interaction cross- 
section, and where 



<$>tot(E v ) 
with 



dt f h dz r 5 ^smt 



L Jo h Jo 



tan6(£,z) 



R 



d + e + z 



■§i ab (E v ,6), 
(7) 

(8) 



B. Large versus Small Ring configurations 

The storage ring geometry is characterized by the length 
of the straight sections D and by its total length L. Be- 
low, we consider the cases of a small (SR) and a large 
(Li?) ring configurations, characterized by (Dsr,Lsr) 
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FIG. 1: Neutrino fluxes scaled by the length of the storage ring L& t ot(E„): The exact results obtained with Eqs.Q-lJBJ with 
a small storage ring SR (solid lines) and a large storage ring LR (long-dashed lines) are shown. The left (right) figure shows 
the fluxes impinging on the small (large) detector (the sizes are given in Tables III and ITU . For the small detector (left), the LR 
result obtained with the analytical estimate Eq. © coincide with the exact result and is not represented here for clarity. For the 
large detector (right), it is also a very good approximation as shown by the dotted- line. The contribution L G &i a b(E v ,9 = 0) 
from the RHS of Eq. © is also presented (dashed lines). All fluxes are obtained with 18 Ne boosted at 7 = 14. 



and {Dlr,Llr) respectively. The results in both config- 
urations can easily be related to one another by split- 
ting the integral over the useful decay path J^ LR = 
f° SR + j£ LR in Eq. Q. Up to trivial l/L factors, the 
LHS corresponds to the LR configuration and the first 
term on the RHS to the SR configuration. The remain- 
ing integral can be given a simple analytical estimate if 
one can neglect the angular dependence of the flux under 
the integral. This happens when the angle under which 
the detector is seen from the extremity of the SR decay 
path ~ R/ (d+ Dsr) is small compared to I/7, i.e. to the 
typical opening angle of the boosted flux. In that case, 
we obtain, for the total flux (J7|: 



^sr 

J LR 



,{$l*{E v )+G$> lab {E v ,6 = Q)}, (9) 



where the geometrical factor G is given by: 



G 



R 2 



AL SR {d + D SR ) 



d + D 



SR 



d + D LR J 



(10) 



The overall factor Lsr/Lrr in © simply accounts for 
the fact that the number of decaying ions per unit length 
is smaller in a larger storage ring, and the second term 
in brackets on the RHS represents the contribution from 
the longer useful straight section. Figure n snows a com- 
parison between the exact flux obtained with Eqs.Q-0 
in both the SR and the LR configurations, and the ana- 
lytic estimate Eq. |(?JJ, for the two possible detector sizes 
considered in the following. We see that the analytical 
formula works very well in the cases considered here. 
Besides, Figure Q] shows that the contribution from the 
longer decay path only brings a ~ 10 % difference for 
the small detector and contributes a factor ~ 2 for the 



larger detector. This already shows that the main dif- 
ference between the LR and SR fluxes comes from the 
geometrical factor Lsr/Llr — 1/15. 

Using the approximate formula for the total fluxes, we 
obtain an approximate relation between the total number 
of events in the LR and SR configurations: 



dN e . 



dt 



L 



SR 



LR 



^LR 



\ dt 



SR 



+ 7 2 (1+ PfGgnh {a}- 



where (<t) 7 denotes the flux- averaged cross section in the 



forward direction 9 = 0: 

r °° dE v $> lab {E v ,e 



0M£„) 



Jo dE v $i ab (E v ,6 = 0) 
Using Eq. the latter can be re- written as: 



(er) 7 = 



r °° dE v <S> cm {E v ) (7( 7 (1 + p)E v ) 
f™dE u $ cm (E v ) 



(12) 



(13) 



It is to be noted that, when the detector is placed close to 
the storage ring, as it is the case here, the total rate © 
depends non-trivially on the geometry of the latter. For 
instance, as discussed above, we observe an approximate 
l/L scaling at fixed D/L in the small detector case. This 
is in contrast with the case of a far detector considered 
in the high energy beta-beam scenarios [23, |32, IH, |34( , 
where the rate is simply proportional to the ratio D/L of 
the straight section over the total length of the ring |44| . 



III. RESULTS 

Here, we present charged-current neutrino interaction 
rates with various target nuclei as obtained from Eqs. 
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FIG. 2: Neutrino fluxes <&tot{E v ) as a function of energy for 
18 Ne nuclei boosted at 7 = 7 and 7 = 14. This corresponds 
to the small ring and small detector configuration. 



FIG. 3: The total rate for the reaction v + D as a function 
of the boost factor 7. This corresponds to the small ring and 
small detector configuration. 



©-© ( Tables HI and ITU) . Four possible nuclei are taken 
as typical examples, namely deuteron, oxygen, iron and 
lead. A detailed study for the case of lead is also done in 
[38J. The "small ring" we consider has 150 meter straight 
sections and 450 meter total length, while the "large ring" 
has 2.5 km straight sections and 7 km total length. The 
detectors are located at a distance 10 meters from the 
storage ring, to allow a maximum shielding of the in- 
duced background in the ring |4£|. For the detector size 
we inspire ourselves on the kinds considered for the pro- 
posed ORLAND facility [ME!- The transverse size is 
chosen so as to catch as much as possible of the boosted 
flux, which main contribution is concentrated in an open- 
ing angle ~ I/7. More precisely, we choose as typical 
dimensions (R = radius, h = depth): R = 1.5 m and 
h = 4.5 m. We also consider the case of a large (kiloton- 
type) water detector with R = 4.5 m and h — 15 m. For 
all detectors here we assume a 100 % efficiency. Finally, 
we have to specify the number of parent ions g injected 
per unit time in the storage ring. According to the feasi- 
bility study HI, 2 x 10 13 6 He/sec and 8 x 10 11 18 Ne/sec 
could be produced with an ISOLDE technique, giving 
about Qd = 10 13 P/sec and g v = 5 x 10 11 v/sec respec- 
tively 30] . An important feature of beta-beams is that 
the number and average energy of neutrinos entering the 
detector depend on the boost factor 7 of the parent ion, 
which can be varied. We present results for two different 
values, namely 7 = 7 (Table I) and 7 = 14 (Table II). The 
corresponding neutrino fluxes are presented in Figure |3 
and range up to about 50 and 100 MeV respectively. 

Let us discuss the number of events shown in Tables 
HJ and UJ The differences between the ^-induced versus 
i/-induced reactions is a combined effect of the relative 
intensities g u /gv = 1/20 and of the different interaction 
cross-sections: the ratio a(y + D)/a(p + D) is roughly 
2 in the whole energy range considered here |13|: from 
UJ, one can see that a(y + 16 0)/a(D + 16 0) is about 



0.5 on average in the energy range relevant to the case 
7 = 7, namely 20 MeV < E v < 40 MeV, and about 
1.5 on average in the range 40 MeV < E v < 80 MeV, 
relevant for the case 7 = 14. The very low rates ob- 
tained for oxygen with 7 = 7 despite the large detector 
size are due to the 15 MeV threshold in the interaction 
cross-section. Next, we observe that the suppression of 
the rates in the LR configuration as compared to the 
SR case for a given 7 roughly corresponds to the geo- 
metrical factor Lsr/ Llr, as expected from the previous 
discussion. In fact, the difference between the LR and 
SR rates can be fully understood by means of the ap- 
proximate relation Eq. tUj . This formula can be used 
to rescale our results for other possible dimensions of the 
storage ring. To this aim, we give the relevant values of 
(er) 7 in each case. When going from 7 = 7 to 7 = 14, the 
neutrino fluxes become more collimated and the typical 
energy of the neutrinos increases. This, together with the 
fact that the neutrino-nucleus interaction cross sections 
rapidly rise with the impinging neutrino energy, increases 
the number of events by more than an order of magni- 
tude. Figure |3 illustrates the rapid rise of the total rates 
with increasing 7. Note that, in the present case, where 
the detector is relatively close to the storage ring, the 
total rates do not have a simple scaling with the detec- 
tor size, due to the non-trivial angular dependence of the 
impinging neutrino flux. 

It is important to emphasize the complementarity be- 
tween low-energy beta-beams and conventional neutrino 
facilities [2(j. The latter provide intense sources of elec- 
tron and muon neutrinos and cover the very low energy 
region, similar to the case 7 = 7 for the beta-beam. Let 
us mention that for comparable neutrino intensities, the 
rates presented in Table [I] are comparable to those ob- 
tained with conventional schemes with detectors located 
at about 50 meters from the source. Low-energy beta- 
beams would produce pure electron neutrino beams and, 
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Reaction 


Ref. 


Mass 


(cr) 7 


Small Ring 


Large Ring 


v+D 


[13] 


35 


36.30 


194 


14 


v+V 


[13] 


35 


23.16 


2494 


178 




[4JJ 


952 


3.33 


60 


6 


^+ 16 o 


[41] 


952 


5.04 


2125 


192 


v+ 56 Fe 


[42] 


250 


137.86 


872 


63 


z,+ 208 Pb 


[17] 


360 


2931.24 


7598 


545 



TABLE I: Number of events per year for 7 = 7 in the small (Lsr = 450 m, Dsr — 150 m) and large [Llr = 7 km, 
Dlr = 2.5 km) ring configurations. These results are obtained by using the exact formulas of Eqs.(6-|BJl. The detector is 
located at d = 10 m away from the ring and has dimensions R — 1.5 m and h — 4.5 m for the D (D2O), 56 Fe and 208 Pb, and 
R — 4.5 m and h — 15 m for the case of le O (H2O), where R is the radius and h is the depth of the detector. The corresponding 
masses are given in tons. The results in the large ring configuration can be precisely understood from those in the small ring 
configuration by means of the analytical formula Eq. 1111 . We give the flux-averaged cross section in the forward direction (a)-, 
see Eqs. G3-G3) in units of 10~ 42 cm . The latter can be used to rescale the present rates for different sizes of the storage 
ring using Eq. I1H . The relevant cross-sections are taken from the indicated references. The results are obtained with 1 year 
= 3.2 x 10 7 s. 



Reaction 


Ref. 


Mass 


(cr) 7 


Small Ring 


Large Ring 


v+D 


[13] 


35 


184.47 


2363 


180 


v+D 


[13] 


35 


96.03 


25779 


1956 


v+ 16 


[41J 


952 


174.28 


6054 


734 


u+ 1(i O 


[41j 


952 


102.00 


82645 


9453 


v+ 56 Fe 


[42] 


250 


1402.11 


20768 


1611 


^+ 208 Pb 


[12 


360 


16310.16 


103707 


7922 



TABLE II: Same as Tabled for 7 = 14. 



by varying the boost factor 7, would offer a unique op- 
portunity to study neutrino-nucleus interactions over a 
wide range of energies. 

To conclude, the present study demonstrates that, 
with typical parameters available from existing studies 
|3fj|. significant interaction rates can be achieved at a 
low-energy beta-beam facility. A small ring - with as 
long as possible straight sections - is the preferred con- 
figuration in the case of a close detector. The rates raise 



rapidly with increasing 7. We think our results are en- 
couraging and we hope they will trigger further investi- 
gations, including, in particular, detailed simulations of 
the detectors. 



We thank J. Bouchez and M. Magistris for useful dis- 
cussions, R. Lombard and M. Mezzetto for careful read- 
ing of the manuscript. 
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[44] For a distant detector (d S> L,D,h), one has simply: 
§tot{Ev) ~ § lab {E v ,9 = 0) #4^, where S = irR 2 is 
the transverse area of the detector. Similarly, one obtains, 
for the rate: ^ ~ g § 7 2 (1 + f3) 2 (a) 7 , where 

Ntarget = n-KB?h is the total number of target nuclei. 



